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Why is natural supersymmetry neither detected nor ruled-out to date? To answer this question 
we use the Bayesian approach where the emphasis in finding prior-independent features within 
broader and minimally biased frames is taken as the guiding principle. The 20-parameter minimal 
supersymmetric standard model (MSSM) global fits to subjective naturalness indicate the existence 
of a prior-independent upper bound on the pseudoscalar Higgs boson mass toa as a function of 
tan/3, the ratio of the vacuum expectation values of MSSM Higgs doublets. For a 30-parameter 
MSSM this implies that rriA < 3 TeV and tan /? < 25 at 95% Bayesian confidence. Removing 
the contradictory subjectiveness within the electroweak fine-tuning measure leads to finding the 
naturalness line, rriA ~ -j= mz tan (3, that reduces by one the number of MSSM Higgs sector free 
parameters. 
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Introduction: Supersymmetry [1] model constructions and phenomenological studies go decades back [2-9] but yet 
have not been discovered nor ruled out by high-energy physics experiments. The specific prediction from supersym¬ 
metry is that there must be new, beyond the standard model, particles with and without colour charges. But it does 
not specify what the particular masses and couplings of the new particles will be. These remain arbitrary with more 
than 100 free parameters. Currently, on the experiments side, it is expected that the large hadron collider (LHC) will 
give a definite answer as to whether low-energy supersymmetry has any role in stabilising the Higgs boson mass at 
125 GeV [10-12]. It is going to probe the so-called “natural” supersymmetry. In view of this we ask the question: Is 
there any robust prediction from natural supersymmetry that could be targeted by the experiments for a discovery 
or an absolute exclusion? We are going to address this question by employing Bayesian statistical techniques. The 
Bayesian method can be used to extract robust predictions from a model based on experimental data. Here the model 
under consideration will be the R-parity conserving minimal supersymmetric standard model (MSSM). 

Robust predictions can be extracted because the Bayesian approach allows a check for the stability of results 
with respect to widely different but well-motivated changes in the prior assumptions concerning the base supersym¬ 
metry parameters. The results which remain the same under the change of the prior assumptions are said to be 
prior-independent and represent the predictions by the model based on the experimental data. Before the LHC com¬ 
missioning, the 20-parameter MSSM fits [13, 14] to neutralino cold dark matter (CDM) relic density, electroweak and 
B-physics data revealed two observables to be prior-independent, namely the (then undiscovered) Higgs boson and 
the lightest top-squark masses. The Higgs boson mass was predicted to lie between 119 and 128 GeV within the 95% 
Bayesian credibility region, while the top-squark mass to be around 2 TeV. The 20-parameter MSSM is specified by 
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where the gaugino mass parameters Mi, M2 and M3 were allowed in -4 to 4 TeV range. The sfermion / mass param¬ 
eters toj vary between 100 GeV to 4 TeV. The trilinear scalar couplings At,b,T,p= e 6 [—8,8] TeV. The Higgs-sector 
parameters rn 2 Hu , rr, /i d , were varied according to m 2 £ sign(m) [—4,4] 2 TeV 2 . The ratio of the vacuum expectation 
values tan/3 = ( H u ) / ( Hd ) is allowed to be between 2 and 60, while sign(ii) the sign of the Higgs doublets mixing pa¬ 
rameter, is allowed to be randomly ±1. The remaining five standard model parameters were also varied in a Gaussian 
manner with central values and deviations according to experimental results [15]. 

In this article we are going to show that by imposing fine-tuning cuts within the 20-parameters MSSM, an additional 
prior-independent result manifests. From this, an inequality relation between the pseudoscalar Higgs boson mass tua 
and tan/3 can be deduced. For the cuts we use the electroweak fine-tuning measure [16, 17] A ew defined as follows. 
Consider the electroweak symmetry breaking condition for a 1-loop corrected Higgs potential, V + AV 
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Here E“ and arise from the 1-loop radiative corrections. For naturalness, each term in the right hand side of Eq. 2 
should be comparable to m 2 z /2 so that 


A ew = maxi ( Ci ) /(m|/ 2 ) 


( 3 ) 


accommodates the fact that for obtaining a natural value of mz then the terms Ci, with i = Hd, H u , p, E“(fc), 
E^(fc), where k denotes the various particles and sparticles contributions, must be of order ro|/2. Using the terms 
that couple the most to the Higgs sector (the case k = ti^, 61 , 2 ) we have 
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The expressions for E^’“(i) are shown in the Appendix. 

In the next section, we describe the Bayesian approach to MSSM naturalness, the fitting procedure and the prior- 
independent result obtained. After that we explain the impact of the result which is a prior-independent bound on 
tua as a function of tan /3 on the a 30-parameters MSSM posterior distribution. We then assess to what extent has 
some relevant 8 TeV LHC supersymmetry limits probe the natural MSSM-30. At the end, we present an analytical 
argument that exposes a subtle methodological contradiction by looking closer at the electroweak fine-tuning measure. 
Fixing the contradiction lead to a no fine-tuning “naturalness line”. After this we summarise our results and give an 
outlook for future studies. 

Naturalness, the Bayesian approaches: There are two major trends in the literature concerning Bayesian approach 
to MSSM naturalness. First, for addressing MSSM naturalness one can compute the amount of fine-tuning at each 
point during the parameters sampling and then penalise highly fine-tuned points according to a chosen subjective limit 
(see e.g. [18]). Within this method, various groups use different fine-tuning measures, e.g. [19-22], The difference 
measures, however, agree when used appropriately as explained in [23, 24], According to the second trend, fine-tuning 
measures manifest implicitly within the Bayesian global fit procedures. In [25-27] it is shown that fitting the MSSM 
parameters in a Bayesian way automatically incorporate a fine-tuning penalisation. Our approach in this article goes 
along the first trend. We use the electroweak fine-tuning measure Eq. 3 and penalise or rule-out MSSM points with 
A ew > 4. The choice A^iy > 4 in search for prior-independent results from global fits to MSSM represents the 
“naturalness” data. A natural MSSM point should have A ew —> 1- Relaxing away from A ew = 1 as a fine-tuning 
cut we choose A ew <2 + 2 where the first “2” represents a 50% fine-tuning and the second a 100% “theoretical” 
allowance on the first. The Bayesian global fit procedure with A ew < 4 is described as follows. 

Fitting procedure: Based on the methodology for our MSSM programme [13, 14, 28-34] the Bayesian global fit of 
the 20-parameters MSSM plus 5 standard model parameters (MSSM-25) were performed separately with linear and 
logarithmic prior probability distributions on the parameters Eq. 1. These were fit to the Higgs boson mass, naturalness 
requirement, neutralino CDM relic density, electroweak and B-physics data shown in Tab. I. MultiNest [49, 50] 
package which implements nested sampling algorithm [51] for exploring model parameters space were used. At each 
MSSM-25 point the supersymmetry spectra were computed via SOFTSUSY [52] and the list of observables O,;, 
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via the following packages. micrOMEGA s [53] was used for computing neutralino CDM relic density QcDMh 2 and 
the anomalous magnetic moment of the muon Sa^-, and SuperIso [54] for predicting BR(B S —>• p + p~), BR(B sq) 
and the isospin asymmetry, Ao _, in B —> K* 7 . With SUSyPOPE [55, 56] we computed the W-boson mass mw, 
the effective leptonic mixing angle variable sin 2 <e;,, the total Z-boson decay width, Tz, and the other electroweak 
observables. These allow the computation of the posterior probability via Bayes’ theorem, 
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Observable Constraint 

Observable 

Constraint 

m w [GeV] 80.399 ± 0.027 [35] 

A 1 = A e 

0.1513 + 0.0021 [36] 

Fz [GeV] 

2.4952 ± 0.0025 [36] 

A b 

0.923 ± 0.020 [36] 

sin 2 9[j f 

0.2324 ± 0.0012 [36] 

A c 

0.670 ± 0.027 [36] 

5a M 

(30.2 ± 9.0) x 10 10 [37, 38] 

Br(B -a X s7 ) 

(3.55 ± 0.42) x 10 4 [39] 

A? 

20.767 ± 0.025 [36] 

Br(B s -a y + y~) 

x 10 -9 [40] 

R° b 

0.21629 ± 0.00066 [36] 

Ra Mb s 

0.85 ± 0.11[41] 

R° c 

0.1721 + 0.0030 [36] 


1.26 ± 0.41 [42-44] 

Afb 

0.0992 ± 0.0016 [36] 

Ao- 

0.0375 ± 0.0289[45] 

Afb 

0.0707 ± 0.035 [36] 

FIcdmIi 2 

0.11 + 0.02 [46] 



m h 

125.6 + 3.0 [GeV] [47, 48] 


TABLE I: Summary for the central values and errors for the Higgs boson mass, the electroweak physics observables, B-physics 
observables and cold dark matter relic density constraints. 
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Here i run over the different experimental observables (data) other than the CDM relic density, x represents the 
predicted value of the neutralino CDM relic density, y = 0.11 is the WMAP central value quoted in Tab. I and 
s = 0.02 the inflated error. The likelihood contribution coming from the CDM relic density is given by Lqdm{x ) 
which is purely Gaussian when the predicted relic density x is greater than the experimental central value y = 0.11 
thus imposing penalisation for CDM over-production. No penalisation is imposed when x < y. The set of experimental 
data used for the fits is 


d. ~ d Tab . I + ew ~ {/+, a i) + {^ew > 50%}. 
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Here d Tab j is the set of experimental central values yt and error shown in Tab. I. T~L in Eq. 6 represents the context 
or hypothesis for the Bayesian theorem, i.e. nature is supersymmetric and that neutralinos make part of the cold dark 
matter relics. From the posterior of the global fits we only show the result which is approximately prior-independent. 
This happens to be an MSSM-25 feature in the (m,A, tan/3) plane. 

Result: The two-dimensional posterior distributions in Fig. 1 shows that requiring fine-tuning no worse than 50% 
as naturalness data while fitting the MSSM-25 to data has a prior-independent impact in the (rriA- tan /?) plane. The 
empty triangular regions are excluded by this naturalness requirement. The prior-independent result is 

m A < ^ tan/3 TeV. (9) 

Eq. 9 is robust and can be applied to any supersymmetry model with a necessary electroweak symmetry breaking 
condition Eq. 2. Next, we assess the impact of this on the posterior sample of an MSSM-30 which favours low values 
of tan /3 < 27 within 95% Bayesian credibility. First we give a brief introduction of the MSSM-30 frame and then 
afterwards check the natural (Eq. 9-based) MSSM-30 points against some LHC supersymmetry limits. 

Naturalness constraint on MSSM-30: In [57], the 30-parameters MSSM was constructed by reducing the parent 
100+ MSSM parameters using a systematic treatment of minimal flavour violation - unlike as done by hand for the 
MSSM-25 case. The parameters consist of e'+Mi, e^ 2 /^, and M 3 in the gaugino sector with Mq 2 (and also their 
imaginary parts which are varied between -4 to 4 TeV. M 3 is allowed to be between 100 GeV to 4 TeV. 

Within the Higgs sector, tua is varied between 100 GeV to 4 TeV while y and Im{y) were allowed within -4 to 4 TeV. 
As for MSSM-25, tan/3 is allowed to be between 2 and 60. The scalar mass and trilinear coupling parameters are 
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The bases Xi,...,i6 are products amongst Kronecker delta S and the Cabibbo-Kobayashi-Maskawa mixing matrix V 
elements. The parameters 01 , 2 , 3 , 6,7 > 0 and 11 , 2 , 2 / 1 , 3 , 6,7 were varied within (100 GeV) 2 to (4 TeV) 2 and — (4 TeV) 2 
to (4 TeV ) 2 respectively; while 04 , 5 , 8 , /m(a 4 , 5 , 8 ), 2 / 4,5 and irn(y 4 , 5 ) were allowed between —8 TeV to 8 TeV. The SM 
parameters are fixed according to experimental results as: mass of the Z-boson, mz = 91.2 GeV, top quark mass, 
m t = 173.2 GeV, bottom quark mass, mb = 4.2 GeV, the electromagnetic coupling, a ~^ = 127.9, and the strong 
interaction coupling, a s = 0.119. The parameters are 

0 ={Mi, 2 , 3 , /x, m A , tan/ 3 , /x), ai, 2 ,...,8,-Bn(a 4 ,5,8), *1,2, 2/1,3,4, 5,6, 7 , 7 m( 2 / 4 ,5) }■ (H) 

The MSSM-30 fits to the Higgs boson mass, the electroweak physics, B-physics, lepton dipole moments and the 
cold dark matter relic density observables disfavour large tan /3 > 30. The corresponding posterior distribution on 
(m A , tan j3) plane is show in Fig. 2(a). The (mA,tan/3) plane is chosen because we aim at showing the impact of the 
prior-independent result Eq. 9 on the MSSM-30 posterior sample. Fig. 2(b) shows what remains after imposing the 
prior-independent naturalness condition Eq. 9 by ruling out the unnatural points. From the surviving posterior, it is 
deduced that m A < 3 TeV and tan/3 < 25 at 95% Bayesian credibility. [71] 

Collider limits on natural MSSM-30 points: To what extent does 8 TeV LHC probe the natural MSSM-30 posterior 
point based on Fig. 9? The natural points can be checked against some LHC limits. ATLAS and CMS 95% confidence 
level limits can be used to constrain models that predict the processes they searched for. Limits on fiducial cross 
sections usually call for writing Rivet [58] analyses to pass over Herwig+-1- [59] Monte Carlo generated supersymmetry 
events. We did not intend to use the full set of such LHC results. Rather, a selected few which are relevant for probing 
the prior-independent naturalness condition Eq. 9 were considered. In [60] a search for scalar particles decaying via 
narrow resonances into two photons is performed. The limits applied on the MSSM-30 pseudoscalar Higgs production 
cross section times branching ratio into two photons did not significantly constrain the posterior sample. The ATLAS 
[61] and CMS [62] limits from search for MSSM Higgs bosons (here the pseudoscalar Higgs) decaying into tau-lepton 
pairs were also considered. These are put next to the production cross section times branching fraction of the 
pseudoscalar decay into tau-leptons for the MSSM-30 posterior computed using FeynHiggs[63]. Fig. 2(d) shows the 
similar case for the ATLAS search for a CP-odd Higgs boson decaying to the Z-boson and the SM Higgs boson which 
in turn decays to tau-leptons [64], All these searches hardly constrain the natural MSSM-30 posterior mostly due 
to the low production cross-section and decay rates of the pseudoscalar Higgs at the LHC. Perhaps, searches with 
topologies involving the pseudoscalar MSSM Higgs boson decaying via charginos and neutralinos could probe better 
the naturalness allowed MSSM-30 region. 

Naturalness line: Here we give a closer look at the numerical and prior-independent result Eq. 9. A zeroth-order 
explanation for the bound on m A as function of tan /3 can be explained using the electroweak fine-tuning measure 
A ew [16, 65]. Consider the electroweak symmetry breaking condition, assuming tan/3 >> 1 but without lost of 




FIG. 1: Two-dimensional posterior distributions from the MSSM-25 fits to experimental plus “naturalness” data. The left-side 
(right-side) plot is for logarithmic (flat) prior fit. The empty triangular region on the plots are explicitly excluded by the 
naturalness limit New < 4. The dashed fine represents the shift that will occur when a relaxed naturalness cut A ew < 20 is 
imposed. The solid contour lines enclose the dark blue (dark) and light violet (light grey to white) regions which correspond 
respectively to the 68% and 95% Bayesian credibility. For both panels dark blue (dark) regions have higher probability compared 
to the light blue to light violate (grey to white) ones. 
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FIG. 2: (a) The posterior distribution for the MSSM-30 on (m A , tan/3) plane, (b) Explicitly shows the effect of the naturalness 
cut Eq. 9 on the MSSM-30 posterior sample. Note that small patches of regions that show up in (b) but absent in the parent 
plot (a) is an artifact of colouring and contouring interpolations/normalisation, (c) and (d) show the effect of the ATLAS and 
CMS 95% upper bounds on the production cross-section times decay branching ratios for the MSSM-30 pseudoscalar Higgs 
boson. For all the panels dark blue (dark) regions have higher probability compared to the light blue to light violate (grey to 
white) ones. 
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Requiring there be no fine-tuning will need all three terms on the right hand side to be comparable amongst 

m 2 

themselves and of order m 2 z / 2. As such tan ^ ~ m 2 z /2 and —m 2 H — [i 2 ~ 0 implies that rn 2 H ^ ~ — fi 2 . In addition 
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( tan^ )/(~ m g„) ~ 1 implies rn 2 Hd » — m 2 Hu . Now applying m 2 Hu ~ —/r 2 and fn 2 Hd » —to the the tree-level 
relation m]\ = 2|/r| 2 + m 2 Hd + gives m 2 A ~ m^. Therefore requiring fine-tuning A ew no worst that A max 


A ew = 


2 rn 2 Hd 
m 2 z tan 2 /? 


< A. 


m A < m z tan/3( A max /2 ) 1/2 . 


(13) 


Loop corrections to the tree-level relation for m 2 A is not going spoil the bound Eq. 13 or Eq. 9. This is the case for 
the loop-corrected [66] m A used for the MSSM-25 fits as can be seen in Fig. 1. There is also no conflict with other 
fine-tuning measures [19-22] since all the measures agree with one another whenever appropriately applied [23, 24], 
As such the core message of this letter goes as follows. We seek for robust predictions for assessing low-energy 
supersymmetry as the model responsible for the Higgs boson mass stability. But this is not possible as long as the 
subjectiveness inherent in the fine-tuning measure Eq. 13 remains. Constructing the bound in Eq. 13 is based on the 
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no fine-tuning and comparability requirements for the terms in Eq. 12. Now imposing A ew A A max is a contradiction 
since this allows fine-tuning even if not worse than 1/A max . Fig. 1 give further insight to this. The fits done with 
A max = 20 shows the corresponding no-go regions similar to the case with A maoo = 4 which do not agree with Eq. 13. 
Our take is that a model point is either fine-tuned, meaning A ew > 1 or not fine-tuned when A ew = 1- This way the 
subjectiveness in selecting a cut on fine-tuning is completely removed. The out come of this is a robust “naturalness 
line” 

mA ~ — -j= mz tan /3. (14) 

V 2 

In fact imposing Eq. 14 reduces the (?7iA,tan/3) plane into a line, meaning one less parameter in the Higgs sector. 
Note that this result is not equivalent with the purely tree-level no fine-tuning measure fi ~ mz■ The ansatz is 
that the naturalness line holds at all loop levels such that radiative corrections to the masses do not spoil the relation. 
The naturalness line Eq. 14 can be used for mapping natural regions of any MSSM frame. 

Conclusions and outlook: We have addressed a question about finding an objective determinant for the existence 
of natural supersymmetry. Our Bayesian approach is based on finding prior-independent features within broader and 
minimally biased frames as the guiding principle [67]. The results of this article and an outlook are summarised as 
follows. 


• The 20-parameter MSSM fits to subjective naturalness, using the electroweak fine-tuning measure, indicate the 
existence of a prior-independent upper bound on the pseudoscalar Higgs boson mass mA as a function of tan fj. 
Imposing the bound on the posterior sample of a 30-parameter MSSM fit to data shows that mA < 3 TeV 
and tan /3 < 25 at 95% Bayesian credibility region. The natural MSSM-30 points are not yet ruled out by the 
8 TeV LHC limits we considered. Constraints from search topologies that include decays into charginos and 
neutralinos could lead to better probe. 

• We seek for robust predictions for assessing low-energy supersymmetry as the model responsible for the Higgs 
boson mass stability. We proposed that this is possible only if the subjectiveness inherent in the electroweak 
fine-tuning measure is removed. Imposing A ew A A max is a contradiction since this allows fine-tuning even if 
not worse than l/A max . A robust method should require their either be fine-tuning, meaning Aew > 1 or no 
fine-tuning, i.e. A ew = 1- This way the subjectiveness in selecting a cut on fine-tuning is completely removed 
and no fine-tuning means mA ~ mz tan/3. We call this relation the “naturalness line”. 

• “Why is supersymmetry not yet discovered?” Up to the public LHC results as of the time of writing this article, 
we claim that an answer is that the regions where it is expected are not yet probed. “Where to look for natural 
supersymmetry?” The proposed regions where it should be expected were derived via Bayesian method with 
minimal model framework construction or theoretical prejudice. Natural supersymmetry should be looked for 
along the “naturalness line” mA ~ mz tan /3 together with a 1-2 TeV lightest top-squarks. The 8 TeV LHC 
limits on gluino and lst-2nd generation sparticles are not in conflict with these predictions. 

• An interesting line for further studies will be to assess the impact of the full set of LHC fiducial cross section limits 
on the “naturalness line” in general and within particular phenomenological frames such as the 30-parameters 
MSSM [57] of the 2-parameters hMSSM [68]. 
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Appendix: The expressions for the E^(ii, 2 , £> 1 , 2 ) contributions to Aew For self-sufficiency we explicitly show the 
expressions for £“(£ 1 , 2 , £< 1 , 2 ) according to [69, 70] 
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where A t = (to 2 ^ — mf )/ 2 + M§ cos2/3(| — |a;w)) g| = ( g 2 + g' 2 )/8, xw = sin 2 dw and F(m 2 ) = 
to 2 (log(m 2 /(3 2 ) — l), with Q 2 = m nip 2 . 2 are computed at tree-level. For the bottom-squarks, 
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where Af, = (m? — m 2 J)/2 + M§ cos2/3(| — Icc^y)- 2 are computed at tree-level. 
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